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We discuss the possibilities of testing the recent claims of a relatively large B — !> D' transition 
form factor through the non-leptonic B —¥ D' + -k~ and B~ — ¥ D'°tv~ decay modes. To estimate the 
width of the latter we need the decay constant of the radially excited D-meson that we computed 
in lattice QCD with Nt = 2 dynamical flavors and after taking the continuum limit we found 
f D i = 117(25) MeV. We also provide an update for the values of the D-meson decay constants 
obtained by using maximally twisted mass QCD. 
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I. INTRODUCTION 



Over the past several years there was a grow- 
ing interest in radially excited D-mesons. BaBar 
collaboration has recently isolated a number of or- 
bitally excited _D-mesons, including a few states 
claimed to be the radial excitations [Ij . For our dis- 
cussion particularly interesting is the 25* state, or 
D 1 meson, to which we focus in the following. From 
the angular distribution of the decay of D' —> D*ir 
the authors of ref. [l| showed that the newly ob- 
served state is consistent with the J p — 0~ assign- 
ment, and since its measured mass, 2539(8) MeV, 
turned out to be close to the quark model prediction 
mw = 2580 MeV [2], they suggested the new state 
could be the 25-state or D'. However, its measured 
width, T(D') = 130(18) MeV, appears to be much 
larger than expected in the quark models of ref. 
which is why the question of identification of the 
new state as a radial excitation remained open. 1 
Note however that in quark models the transition 
matrix elements to radially excited states are hard 
to control because the position of the node in the 
wave function of the radial excitation highly de- 
pends on the type of equation used to derive wave 
functions. 

Besides the mass and width of the radial exci- 



tation it is particularly interesting to examine its 
weak interaction properties. Recently, in ref. Q, 
it was suggested that a potentially large B(B — > 
D'tu) could help solving the "1/2 vs. 3/2 puzzle", 
because the subsequent Y(D' — s- Di/27r) is much 
larger than T(D' — > D^^n) due to the fact that the 
emerging pion is in its s- and d-wave respectively. A 
large B(B — > D' Iv) would therefore result in an ex- 
cess of the detected B — > D 1 / 2 {'x)£v with respect to 
B — > Z?3/ 2 (7r)^V. In a simplified notation, by -D1/2 
we labelled a state belonging to the j p = (l/2) + 
doublet of mesons [Dq,D[], while -D3/2 stands for 
a meson from the j p = (3/2)+ doublet [Z?!,!?^]. 
The "1/2 vs. 3/2 puzzle" refers to the fact that the 
B- factories observed B(B — > D 1 / 2 lv) ~ B(B -» 
D^/2^v) which turned out to be in conflict with 
redictions B(B ->■ D 1/2 tv) < B(B -> D 3/2 £v) @- 



For a recent review on the orbitally and radially charmed 
mesons, please see ref. 



An argument in favor of significantly large 
B — >• D 1 transition form factor has also been ad- 
vanced in ref. [10] as it would provide a part of 
the suppression to the B — > D^*> form factors 
that could then explain the tension between |V c f,|, 
Cabibbo-Kobayashi-Maskawa (CKM) coupling, ex- 
tracted from the inclusive and that extracted from 
the exclusive semileptonic decays. 

In this paper we discuss the possibility to ex- 
perimentally test the recent claims of significantly 
large T(B — > D 1 iv) by measuring the corresponding 
non-leptonic decays. That possibility was first men- 
tioned in ref. @ which we further develop here. By 
using the factorization approximation, which ap- 
pears to be very well verified by the data in the 



2 



case of the so called Class-I decays [111], one can 
compare B(B -> Dir) with B(B -> Z?'7r) and de- 
duce 



B(5° 



D+7T 



(1.6 ± 0.1) x |/£(0)| 2 , (1) 



so that for a large value of the B — >• £)' form fac- 
tor, /<f (m 2 ) w /« (0) = /f (0), one would have a 
large number of events that should be detectable in 
the data-sample accumulated at BaBar and Belle, 
and especially in that of LHCb. Furthermore, by 
means of numerical simulations of QCD on the lat- 
tice (with N{ — 2 dynamical light quarks) we com- 
puted the decay constant fjji and obtained 



h 

f. 



= 0.57 ±0.16. 



D 



(2) 



so that the ratio between the Class-Ill and Class-I 
decays becomes 

2 



B(B - -> D'°tt- 
B(B 



D'+ir 



Tb- 



1 



0.13(4) 



(3) 



The above estimates are made by assuming that 
the radial excitation is indeed the state observed 
by BaBar, raw = 2539(8) MeV. As we shall see, 
the value of mw computed on the lattice is larger. 
That however results in only slightly modified num- 
bers on the right hand side (r.h.s.) of eqs. (|l|3j) . and 
our statement that the B — > D'it decays should be 
detectable in the _B-experiments for non-negligible 
/+ (0) remains valid. 

As a by-product of this study we improved 
the estimate of the Z?-meson decay constants 
made by the European Twisted Mass Collabora- 
tion (ETMC) [H, El to 



f Ds = 252 ± 3 MeV, 



Id 



= 1.23(2)(1), (4) 



which translates to f D = 205(5) (2) MeV. 

The remainder of this paper is organized as fol- 
lows: in Sec.|IT]we introduce the form factors, com- 
ment on their existing estimates and describe how 
they could be extracted from the non-leptonic de- 
cays; in Sec. [TTT] we discuss the lattice determina- 
tion of several quantities including fu< and fr>> ', in 
Sec. IIVI we complete the phenomenological discus- 
sion of Sec. Ulland finally summarize in Sec. |V] 



II. HOW LARGE IS THE B -> D' FORM 
FACTOR? 

The weak interaction matrix element driving the 
B — > D' ' iv decay is conveniently expressed in terms 



of hadronic form factors /(f + (g 2 ) as, 



{D'(k)\c^b\B(p)) = 



9 9 

m o — mi 



^ I?(v 2 ) , 



(5) 



where q = p — k, and the superscript "i?" is used to 
distinguish B — > D' transition from the more famil- 
iar B — > D weak decay form factors. At q 2 = the 
two form factors are equal, (0) = f$ (0). Contri- 
bution of the form factor /^ (q 2 ) to the differential 
decay rate of B — > D'lv in the Standard Model 
comes with a factor oc mj, due to helicity suppres- 
sion, and for the case of t = e, fi it can be safely 
neglected. 

Form factors are nonperturbative quantities and 
for a decay to the radial excitation they are very 
hard to extract from numerical simulations of QCD 
on the lattice. Apart from an early attempt in 
ref. [14] . there are no results concerning the form 
factors /(f + (q 2 ) reported so far. Instead, there 
are estimates made by using the constituent quark 
model QjKIll- An attempt to extract /+(0) in the 
framework of QCD sum rules has been reported in 
ref. @. 

Similarly toB^ Dlv decay, also in this case the 
considerations in the heavy quark limit are very 
useful. For m c .b — > oo the form factor I+(q 2 ) is 
related to its Isgur-Wise function £r(w) as, 



mB + mo' 



(6) 



where q 2 and the relative velocity w are related 
via, q 2 = m 2 B + m 2 D , — 2mBT)iD'W. Since the 
ground state and the radial excitation are orthogo- 
nal, the Isgur-Wise function at the zero recoil limit 
is £r(1) = 0. This is in contrast with the usual 
Isgur-Wise function for the elastic transitions [ljj 
in which case £(1) = 1, and £(u>) is then conve- 
niently parametrized by the first few terms in ex- 
pansion around the zero recoil limit, 

t{w) = l-p 2 {w-l)+ l -o 2 {w-l) 2 . (7) 

Stringent bounds on the shape of the Isgur-Wise 
function £(w) were derived in ref. (l8j . of which the 
most interesting one for our purpose is the relation 
between the slope p 2 and curvature tr 2 of the ordi- 
nary Isgur-Wise function on one side, and the slope 
of the Isgur-Wise function relevant to the transition 
matrix element between the ground state and the 



3 



radially excited one £r(w), namely 



1/41 = I^WI <\ll^-\p 2 -{p 2 ) 2 - (8) 



Experimenters fit their B 
form [H, 



Dfiiy data to the 



G(w) 
G(l) 



'D 



z + {51p 2 D - 10)z 2 + . . .(9) 



where z = {y/w + 1 - y/2)/(y/w + 1 + a/2), so that 
p 2 and a 2 in eq. ([7]) can be identified as 



P =PD > ° = 32 



(10) 



which, together with the measured p 2 D 
1.19 (4) (4) [Hi, inserted in eq. © gives 



1/41 < o. 



(11) 



Since £,r{1) — by definition, the only plausible 
way to enhance /+ (g 2 ) is through large power cor- 
rections at the zero recoil point (w — 1), i.e. to 
/+(<7ma X ) wher e Q^ax = ( m B - m D ,) 2 . In this 
way the form factor would be significant in the low 
q 2 region too which due to the large phase space 
enhancement would provide a significant B(B — > 
D'lv). The numerical estimates of this form fac- 
tor have been made in the framework of two quark 
models, and both values are in fact large: /f (0) = 
0.37 - 0.41 [H[, and (0) ~ 0.2 [16]. A rough 
QCD sum rule estimate gives /£ (0) = 0.16(11) [|. 



A. Non-leptonic decays can help 



From the above discussion we see that a large 
value of f+ (q 2 ) at low q 2 's is difficult to reconcile 
with the bound {5} unless huge power corrections 
modify the heavy quark mass limit of the form fac- 
tor. The size of /+ (0) can be tested experimentally 
by considering the Class-I non-leptonic decays for 
which the factorization approximation is known to 
be very well verified in experiments and is for mally 
proven to be valid in the heavy quark limit (22l. |23| . 

More precisely, we will consider B — > D' + n~ de- 



2 Radiative corrections to the bounds of ref. [ill have been 
computed in heavy quark effective theory in ref. [Tgj ] and 
are found to be very small. 



cay for which the factorization amplitude reads, 
Gf 



A 1 



-prO, 



fact 



V2 



V cb V: d a^~\dl^u\Q){ D'+\c^b\B ) 



i— ^V c bV* d aiU(m 2 B 



m 2 D ,)f R (ml), 



(12) 



where di(mb) is the Wilson coefficient that contains 
information about physics at short distances, and 
/rf (m 2 ) « / Q fl (0) = /f (0). When combined with 

the more familiar B — > D + tt~ decay mode, one 
gets 



B(B° 



m 



D' 



D + 7T-) 



m 



D 



X(m Bl m D , m v ) 



1/2 



/+(0) 



(13) 



where X(x, y, z) — [x 2 



(y + z) 2 }[x 2 - (y 
and /+(0) is the B — >• D form factor, the value 
of which has been measured experimentally. In 
the lowest bin, m 2 < q 2 < 1 GeV 2 , the BaBar 
collaboration obtained \V c b\f+(q 2 ) ~ |V^f,|/+(0) = 
0.02642(8) |24j, which then, after using \V cb \ 
0.0 111 ; If.: gives /+(0) = 0.64(2). With this 
last number, and with mw = 2.54 GeV, eq. ([T3"|) 
gives 



6(5° 



(1.65 ±0.13) x |/^(0)| 2 .(14) 



S(S U 



D+TT- 



B° 



D+ir- 



has been extensively studied at the 
S-factories, and the final result B(B — > D + tt~) = 
0.268(13)% HI] can be combined with eq. (0 to 
predict 



6(5° 



\f+(0)\ 2 x (4.7 ±0.4) x 10~ 3 , 

(15) 



which for the values of /+ (0) discussed in the litera- 
ture leads to B(lf -> D'+ir-) - 10~ 4 , that should 
be possible to measure even in the sample of BaBar 
and Belle, and especially at LHCb. 

If the experimenters indeed detect a significant 

number of B —> D' + ir~ events then it would be 
interesting to check if the picture is consistent with 
what happens in the Class-Ill decay B~ — > D'°tt~ . 
The factorized amplitude of this mode contains two 
pieces, 



All I 

•^fact 



1/2 d 



aifn(m 2 B 



m 2 D ,)f R (ml) 



+a 2 f D ,(m 2 B -ml)f^(m 2 D ,) 



(16) 
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where /(P^^^f)/) is the B — > -kIv decay form fac- 
tor. One could then consider, 



B[B~ 



D 



>0o 



-> D'+TT-) 



T B - 

rgo 



1 , ° 2 

H x 

Ol 



mi 



/o 



/f(0) 



/d /tt 



(17) 



in which two unknown quantities are fw / Jd and 
the ratio of the Wilson coefficients a^jai. The first 
quantity can be determined on the lattice, while 
the second one can be extracted from the measured 
B(B~ -> D°ir-) = 0.481(15)% and the above- 
mentioned B(jf ->■ D+TT-) = 0.268(13)% HI]. 



III. LATTICE DETERMINATION OF THE 
DECAY CONSTANTS f D(s) AND f D , 

In this section we present the results of our com- 
putation of masses and decay constants of £)W- and 
Dg -mesons by using the maximally twisted QCD 
(MtmQCD) action on the lattice [23] with N f = 2 
dynamical light flavors. At fixed lattice spacing, we 
will also compare the MtmQCD values with the re- 
sults obtained by using the standard Wilson-Clover 
action with N{ = 2 dynamical quarks and with 
those obtained in quenched QCD. Moreover, we will 
improve the values for /d (s) and /_d (s) / /d reported 
in refs. did. 

Extraction of the mass and the decay constant 
of a given hadron state on the lattice is made from 
the study of the two-point correlation function, 



X 

= -(^Tt\r i S e (0 > x)T J S q (x,0)]), 



(18) 



where Or ; = cTiQ is the bilinear quark operator, 
with c and q being the charm- and the light-quark 
field respectively, and Tij chosen to ensure the cou- 
pling to the state with desired quantum numbers. 
In our study q is either the strange quark, or it 
coincides with the light sea quark. In the above 
notation S q (x, 0) = (q(x)q(O)} is the quark propa- 
gator computed in the background gauge field con- 
figuration by inverting the Wilson-Dirac operator 
of MtmQCD on the lattice. The simplest and the 
most convenient choice of the operators is to use 
Or; j — P5 — C759, and extract the mass and decay 
constant of the lowest lying state from the exponen- 



tial fall-off of (fT5|) which for large time separations, 

t > 0, 



c 55 (t) = C£p 5 (x; 0^(0; o)> 



1 2 cosh[m Dg (T/2 - t)} „_ mDnT/2 

\ ' 

(19) 

with T being the size of the temporal extension of 
the lattice, and Zo q = {§\c^$q\D q ). In cq. (fT9|) 
we used the symmetry of the correlation functions 
with respect to t O T — t of our periodic lattice. 
To extract the radial excitation properties one can 
subtract the r.h.s. of eq. (fl~9f from the correlator 
C 55 (t), 



C' 55 (t) = C 55 (t)- 

, 2 cosh[m Dq (T/2-t)} 



-m D T/2 



(20) 



cosh 




c 55 (t) 


cosh 


mf, q {t) (^-t-l 


) 





and check whether or not there is a plateau of the 
effective mass, m c p, (t), defined as 



, (21) 



and possibly fit to the form similar to eq. (|T9|) to ex- 
tract the mass and the decay constant of D' q . This 
strategy can be extended and combined with corre- 
lation functions computed by using different source 
operators. 

Another way to proceed is to work with several 
interpolating operators that can be easily built if in- 
stead of the local fields q and c one uses the smeared 
ones, q ng and c„ s , that we generically call ip ng , de- 
fined via 



1 



KgH 



1 + 6k„ 



where the smearing operator H reads [26 



^=1 



r/n»T x. 



(22) 



(23) 



with being the n a times APE smeared link [29[ , 
defined in terms of (n a — 1) times smeared link 
t/^™" and its surrounding staples of links denoted 
by K 



Proj SU( 3) [(l-a)U^ + ^V^ 



(24) 
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The above steps are known as the Gaussian smear- 
ing procedure. In this work, we choose the following 
values of the parameters: 

k cj = 4, n g e (0, 2, 10, 32), a = 0.5, n a = 20 . (25) 

We checked that the correlation function computed 
with both quark fields smeared is equal to the one 
obtained with only one field smeared but with twice 
as many smearing steps, n g . However, we observe 
that the correlation functions computed with both 
fields smeared are less noisy and for that reason 
the results presented in this work are obtained by 
using both q and c fields smeared. Therefore, with 
various values of n g , and for T = 75 we get various 
operators Oi that can be combined in the matrix 
of correlation functions (fT8|) . Note that the choice 
n g = corresponds to the local operator which is 
needed for the computation of physically relevant 
decay constants. The problem of extraction of the 
hadron masses from 



C l3 {t) = (O^)Oj(O)} 



7 ■ 7* 



cosh 







-)] 


m (n) 


(§- 









is then reduced to the generalized eigenvalue prob- 
lem (GEVP) m 



C l3 {t)vf\tM) = A(™)(Mo)a;(*oK n) (Mo),(27) 

where Cy(to) is chosen for computational (numer- 
ical) convenience, while \ {n) (t,t Q ) and v^ n) (t,t ) 
are the eigenvalues and eigenfunctions of the ma- 
trix C _1 (to)C(t). The goal we achieve by solv- 
ing GEVP is the identification of interpolating field 



.(») 



Oi, which has a desired property that 



(£>W|O(»)t|0) 
tained from 



A (n) (Mo) = 



A n &mn- Masses are then ob- 



(28) 



cosh 


m (r 


,(T/2- 


■*)' 


cosh 


m in 
iJ <l 


(T/2- 


to) 



where n = 1 corresponds to the lowest lying pseu- 
doscalar D q mesons, and n = 2 to their first radial 
excitation D' q . To extract the decay constant one 
needs the matrix element of the local operator and a 
state \D q n ^) isolated by using O'™'. In other words, 



(D^\O{\0) = 



i 



(29) 



In the case of MtmQCD on the lattice, the local 
operator of interest is Ol = P5 = cj^q because 
(fj,q + fi c )P5 is renormalization group invariant, and 
therefore no renormalization constant is needed to 
compute the pseudoscalar decay constant. 3 This is 
not so in the case of Wilson-Clover action where it is 
more convenient to use Ol = Aq = -^(So) •270759) 
with ^((7g), the axial current renormalization con- 
stant. However, since we shall be interested in the 
ratio of the decay constants, fr>' /fo q , one can use 
Ol = P5 in the case with the Wilson- Clover action 
as well. 



A. Lattices used in this work 

The lattice results of this paper are obtained by 
relying on the ensembles of gauge field configura- 
tions produced by the ETMC [3l| from the simula- 
tions of MtmQCD [27j with Nf — 2 dynamical mass 
degenerate light quarks. Main information concern- 
ing these ensembles is given in Tab. HI We use the 
results of ref. 32J to fix the charm (c) and strange 
(s) quark masses at each lattice spacing and then 
compute the correlation functions needed for the 



(26) extraction of strange and non-strange D^-meson 



properties. The quark propagators are computed 
by using stochastic sources, and in the computation 
of the correlation functions we used the so called 
one-end trick 31 1. 

The effective mass plots are obtained from 
corresponding eigenvalues \( n \t,to), solutions of 
cq. (|2"7|) . which by virtue of eq. ([2"5f can be writ- 
ten as, 



m^f„) (t) = arccosh 



A< n >(t + l,to) + A< n >(t-l,to) 



2A(»)(Mo) 



(30) 



For the case of the lowest lying state (n — 1) and 
the first radial excitation (n = 2) the signals are 
illustrated in Fig. Q] In the plateau region, each 
m /.At) is then fitted to a constant m„(o. We 

checked that the results for the first radial excita- 
tion remain stable when we change the size of the 
matrix of correlators. We also checked that from 
this study we cannot extract a signal for the second 
radial excitation: only a few points at t < 5 can be 
seen before the error bars become overwhelmingly 
large. We also checked that the choice of to in the 
GEVP (|2~T|) does not make any impact on the re- 
sults presented here. Finally, we also note that the 



fiq.c is the quark mass parameter. 
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FIG. 1: Effective mass plots leading to the masses of D 3 (left column) and D' s states (right column), deduced 
from the matrix of correlation functions as discussed in the text. Plots are provided for all four lattice spacings 
considered in this work with MtmQCD action. 



radial excitations extracted on the plateaux of sub- 
tracted correlation functions ([20)1 are completely 
consistent with those obtained from the solution 
to the GEVP [g7J). 

Concerning the decay constants they are ex- 
tracted from the matrix element obtained by using 
eq. (|29|) . with Ol = P5, and the definition 



\P^\Df)=m 2 D „ ) f D „ ) . CM) 



The fitting intervals for extracting the masses and 
decay constants for the lowest lying states are 

t/ae [8, 22^=3.8, t/ae [8,22^=3.9, 
t/a £ [12, 26]^ =4 .05, t/a e [14, 30)0=4.2, 

while for the radially excited states the following fit 
intervals have been chosen, 



t/a € [6, 10^=3.8, 
t/a e [9, 12^=4.05, 



t/a e [9, 12^=3.9, 
t/ae [11,14^=4.2. 



For some values of the sea quark mass we have a 
few more points to fit but globally the time intervals 
noted above are used to obtain the results that we 
present in tabs. Illll and HVi in lattice units and for 
each of the lattice setups employed in this work. 



B. Re-evaluation of fo s and /d s //c 

The results of ref. |l2[ included the simulations 
at three different lattice spacings and the value 
f Ds = 244(8) MeV has been reported. That value 
has been improved in ref. [l3j where the simulations 
at a smaller lattice spacing have been included in 
the analysis, leading to f Ds = 248(6) MeV. Fur- 
thermore, while improving the MtmQCD estimate 
of fD s /fD, the authors of ref. [H[ also added the 
systematic uncertainty related to the chiral extrap- 
olation, which was omitted in ref. 12]. Their fi- 
nal result, JdJ/d — 1-17(5), allowed to deduce 
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f D = 212(8) MeV. 

Results in refs. [H], [l3| have been obtained from 
the correlation functions with local source opera- 
tors only. In the present work we implement several 
levels of the smearing procedure discussed above, 
with parameters (|25|) , and then combine the result- 
ing correlators in a matrix. Solution to the GEVP, 
together with a slightly modified procedure to ex- 
tract fu s and fuj Id, result in more accurate re- 
sults which is why in this subsection we update the 
values presented in refs. [l2l Il3j. 

To get the physically relevant fp s we need to 
extrapolate the values we obtained from all of our 
lattice ensembles (c.f. Tab. lIVp . We choose to com- 
bine fo s and rn,D s in the dimensionless ratio that 
we then fit to the form, 



h 



niD s 



latt. 



A, 



1 + B D m q + C Ds 



0/3=3.9 



(32) 



where Ad„, Bd 3 and Cd s are obtained from the 
fit and their values are given in Tab. [Vj The above 
form takes into account the fact that the lattice dis- 
cretization effects to the hadronic quantities com- 
puted in MtmQCD are cx a 2 . Note that we di- 
vided by the lattice spacing obtained at (i = 3.9 
so that the parameter Cd s actually indicates a size 
of discretization effects at j3 = 3.9. After taking 
m q = mf s (2 GeV), also listed in Tab-HTD the fit of 
our data to eq. ([32]) in the continuum limit and at 
the physical mjff (2 GeV) = 3.6(2) MeV *32j, gives 



D. 



mo. 



ph. 



0.1281(11) 



(33) 



With the help of rnP D ' 
finally have 



1968.5(3) MeV [H, we 



f Ds = 252(3) MeV. 



(34) 



We checked that this result remains stable if we 
omit from the continuum extrapolation the results 
obtained at /3 — 3.8. This result is also consistent 
with those obtained from simulations with Nf = 
2 + 1 flavors of staggered quarks in the continuum 
limit [33j, with those computed with N{ = 2 + 1 
flavors of Wilson- Clover quarks at the single lattice 
spacing pjij , as well as with the recent experimental 
results presented in ref. 135?] . 

As for the SU(3) light flavor breaking, the ratio of 
Jd 3 / Jd is combined with the meson masses in the 
way consistent with the heavy quark expansion, 



<>D a 



(35) 



where the index u q" labels the valence light quark, 
which in our study is mass degenerate with the sea 
quark. As in ref. [T^, EH we fit our results to a form 



„latt. 



— Aj 



1 



x 3 - 1 



3 5 5 



4 (4tt/) 5 



ml log(m 2 ) 



+B r mi + C r 



0/3=3.9 



(36) 



where for X — we have the expression similar to 
the one used in eq. (|3"2"|) . and for X = 1 the extrap- 
olation formula includes the chiral logarithmic cor- 
rection that has been computed in the framework 
of heavy meson chiral perturbation theory [36j . To 
use the latter formula one needs to fix the value of 
the soft pion coupling to the doublet of the low- 
est lying heavy-light mesons, g. That coupling has 
been recently computed in ref. [37j on the same sets 
of gauge field configurations that are used here, and 
the result is g = 0.53(3) (3). The results of the fit of 
our data to eq. (|36l) are collected in Tab. [V] Here 
we note that, 

forV = 0, 4^ = 1.128(10), 

3D 



for X = 1, 



fD 



1.227(13). 



(37) 



After averaging the last two results, we finally have 



= 1.177(16)(50), 

Id 



(38) 



from which we can deduce fjj = 214(5) (9) MeV, 
where the second error reflects the systematics 
arising from the chiral extrapolation. In order 
to circumvent the large logarithmic correction in 
eq. (f3"6"l) . one can study a double ratio [38j . 



R a 



mp 3 fp s I fp q 

m Dq fKa/Uaa 



(39) 



for which the logarithmic term is about 10 times 
smaller than in eq. (I36|) . 



R 



latt. 



Aj 



l + X 



9g 2 



4 (4^/) 2 
+B R ml + C R 



ml log(m 2 ) 



09=3.9 



(40) 



and therefore the difference between the values ob- 
tained by setting X — 1 and X = is much smaller, 
which can also be appreciated from the plot shown 
in Fig. [2j We get, 



^ = 0.995(9)(4) x Ip- 

JD U 



(41) 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 
m n 2 [GeV 2 ] 

FIG. 2: Chiral extrapolation of the ratio of decay constants: Upper curves correspond to the fit of our data for 
4>d„/4>d according to eq. i|36l) with X = and X = 1; Lower curves correspond to the chiral extrapolation of the 
double ratio (4>d b / 4>d) / (Jk //it) with the formula i4Lfy with X = and X = 1. The vertical line indicates the 
physical pion mass. Dashed (solid) curves depict the extrapolation without (with) inclusion of the chiral logarithms. 
The symbols corresponding to the lattice data are: O for j3 — 3.80, ■ for j3 — 3.90, • for j3 = 4.05, and □ for 
f3 = 4.20. Note that the result of the linear extrapolation of the double ratio has been slightly off-set to make it 
distinguishable from the one in which the chiral logarithms have been included. 



where the central value is obtained by averaging the 
results of extrapolations with X = and X = 1, 
and the second error reflects the error due to chiral 
extrapolation. The results of the fit of our data 
to eq. (l4Tfl) are listed in Tab. [V] Following the same 
strategy as in ref. [TH, from the results for fn q / }-n qq 
listed in Tab.[TV]we obtain f K /U = 1-23(1), which 
then gives, 

^= 1.23(2)(1), (42) 
Id 

that combined with fr> a in eq. (|34l) gives 

f D =205(5) (2) MeV. (43) 

We note also that the above result remains remark- 
ably stable if the data on our coarser lattices (cor- 
responding to (3 = 3.8) are left out from the chiral 
and continuum extrapolation. 



C. Ratios f D 'JfD s and fo'lfn 

We now discuss the masses and the decay con- 
stants of the radially excited D-mesons. We fo- 
cus to the dimensionless mr}' q /mr) q and fD' q /fD q , 
that are easily built from our results presented in 
Tab. Mil (non-strange) and Tab. IIVI (strange). In 
the following we denote by T one of the four quan- 
tities discussed in this section, namely mo'/mr). 



m D >Jm Ds , Jd'/Id, and /dj//d„, and fit each to 
the form similar to eq. (|32|) . 



-putt. = A 



1 + Bjrrriq + C T 



Oa=3.9 



(44) 



We get the following physically relevant results, 
1.53(7), M = o.53(9), 

IDs 



^ = 1.56(9), f -f 



0.50(12). (45) 



An illustration of that fit in the case of m £>< /to_d 
and fjji I ju is provided in Fig.[3J while the values of 
Ajr, Bp, and Cj= for all four quantities can be found 
in Tab. [V] We observe that the above ratios do not 
exhibit a regular behavior in a 2 , and are practically 
independent of the light quark mass. For that rea- 
son it is tempting to fit our data to a constant, i.e. 
to impose Bjr = Cjr — in eq. PH) . We obtain 



mD ' s 
m D 



1.52(2), 



D' 



ID, 



1.56(9), ^ = 0.65(3). (46) 

JD 

In other words the mass ratios remain stable while 
the ratios of decay constants change quite consid- 
erably. We decide to take the difference between 



0.65(3), 
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2.0 r 



S 



1.5 = EiUlII=^ 




0.0 



10 20 30 40 
m q [MeV] 



50 



60 



FIG. 3: Chiral extrapolation of the ratios mo'/mD and /d'/Id by using eq. Dashed lines corresponds to 

central values of various lattice spacings, and the bands indicate the error bars of extrapolation. The symbols of 
the lattice data points are the same as in Fig. 



the central values in eq. (|45j) and in eq. (|46j) as an 
estimate of systematic uncertainty and after sym- 
metrizing the error bars we finally obtain: 



mD ' s 



1.53(7), ^ = 0.59(11), 



1.55(9), ^-=0.57(16). (47) 

JD 



With respect to the mass of the state suggested 
to be interpreted as D' and measured at BaBar, 
mo'/mo — 1.36, our result is significantly larger. 
If the state measured at BaBar is indeed the radial 
excitation D', the fact that our value is larger would 
be difficult to explain. 

One could suspect that tuning the twisting angle 
to its maximal value on the lattice could be made 
only up to discretization effects which then induce 
a pollution to the extraction of the desired hadron 
state by the state with opposite parity. However, 
since that pollution is an 0(a 2 ) effect [27j and since 
we perform the extrapolation to the continuum 
limit, that argument could not be used to explain 
the potential discrepancy between our findings and 
the value measured at BaBar. 

To further study this issue, we used the data at 
fixed lattice spacing (corresponding to (3 = 4.05) 
and compared them with the results obtained by 
using the Wilson-Clover quark action with Nf = 2 



light flavors at nearly equal lattice spacing. 4 The 
results are 

MtmQCD : — ^ = 1.55(6), ^ = 0.69(5), 
Clover : — ^ = 1.48(7), ^ = 0.77(9). 

m D s JD S 

In other words, at fixed lattice spacing and fixed 
light sea quark mass [mf s (2 GeV) ~ 25 MeV], the 
results for mt3> /rriD a and fo' / fo s obtained by us- 
ing MtmQCD action are consistent with those ob- 
tained with the Wilson-Clover action. We therefore 
conclude that, within the above error bars, the mass 
of the radial excitation obtained on the lattice with 
Nf = 2 dynamical quarks is larger than the state 
measured by BaBar. 

Another potential difficulty when computing the 
properties of higher excited states on the lattice 
with light dynamical quarks is that the decay chan- 
nels with emission of a pion might open up and 
modify the value of extracted mass and decay con- 
stant. In the problem at hand, such channels 
are D' -» D*n and/or D' -> D^ir in the case 
of non-strange radial excitation, and D' s — > D*K 
or D' s — > D* Q K for the strange radial excitation. 
That difficulty does not exist in quenched QCD 



4 The gauge field configurations with Wilson-Clover action 
have been produced by the QCDSF collaboration [39( that 
we have used in our previous work 40]. See also Tab. ITT1 
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(Nf = 0) which is why we produced a set of 
quenched QCD configurations at a similar lattice 
spacing (c.f. Tab. [TTJ) and computed the mass and 
decay constant of D s meson and its radial excita- 
tion by following the same steps as described above. 
We have 



N f = 
N f = 2 



m D ' a 
rriD s 



1.41(9), 



fD. 



1.48(7), S -Bk 

IDs 



= 0.67(12), 
= 0.77(9). 



From our data we cannot see the effects of the sea 
quark mass. We therefore conclude that within 
the statistical errors of this study the radially ex- 
cited D' state is heavier than the one measured 
by BaBar as well as the onepredicted by the con- 
stituent quark model of ref. 0] . That issue deserves 
more research, especially exploring the simulations 
with very light sea quark mass, in order to check 
whether or not the dependence of mrj/ on the light 
sea quark mass changes considerably close to the 
chiral limit, similar to the findings of ref. 41| for 
the Roper resonance. An indication that this in- 
deed could be the case is provided by the results 
reported in ref. [42j . where the simulations at one 
lattice spacing have been used to compute the spec- 
trum of _D( s )-mesons. J 

Concerning our results for the decay constants 
of the radial excitations, we see from eq. (l4"5t that 
they are considerably smaller than the those of the 
lowest states. This situation is qualitatively differ- 
ent from what happens in the heavy quark limit 
(m c — > oo), in which [44[ 



fp' q 



lim 



> 1 



(48) 



while in our case, with the propagating charm 
quark, the above ratio is smaller than one. 



IV. MORE PHENOMENOLOGY 

With the results of our lattice study given in 
eqs. (|43|45|) . we now return to eq. (JTTj) in which 
the only remaining ingredient we need is the ra- 
tio of Wilson coefficients 02/01. Their values have 



4J. 
but 



5 Very recently two preliminary studies appeared in ref. 
They do not report the numerical values for m u i 

from the plots provided we could see that their m D i is 
larger than the mass of the state observed at BaBar, 
2538(8) MeV. 



been computed through matching of the full Stan- 
dard Model with the low energy effective theory in 
which the properly resummed next-to-leading loga- 
rithmic QCD corrections have been included. The 
final result in the MS(NDR) scheme at /i = m& is 
a2/ai= 0.172 45]. In the factorization approxima- 
tion, however, 02/01 is considered to be an effective 
parameter that can be extracted from similar non- 
leptonic decay modes [III US E3| • We can use the 
corresponding B — > Dtt decays, and from 



B(B~ -> £>°7r-) 

B(~lf ->■ D+TT-) 



T B - 
Trc 



1+ - 



eff. 



f B 
JO 



'"(ml) 



/r^(o) 



m| - 


ml 


m| 




m 2 D 




2 




U. 







we get 



— I = 0.368(65) . 



eff. 



(49) 



(50) 



To get the last number we used the following in- 
put values: (i) g(if ^ D+ti-) = 0.286(13)%, 
B(B- -> D°tt-) = 0.481(15)% Hi, (ii) TbJt b - = 
1.079(7) [11, (iii) ff D {0) = 0.64(2) [H, (iv) 
f B ^(m 2 D ) = 0.29(4) 0, (v) f D /f„ = 1.56(3)(2), 
obtained in this work, c.f. eq. (|4T)) . 

Therefore, with our fw/ Id = 0.57(16) inserted 
in eq. (ITT)) , we have 



B(B- 



B(B 



D'+tt- 



tb- 



0.14(4) 
IfW 



(51) 



where we used our result mw /mo = 1.55(9). If, 
instead, we used m^' /mo = 1-36 as suggested by 
BaBar, the numerator on the r.h.s. side would 
change to 0.13(4), as quoted in eq. ([3]). 

In addition, by using our value for mjy/mo the 
result of eq. (TT3"j) would become, 



B(W 



D 



'+0 



B(B° 



D+TT- 



(1.24 ±0.21) x |/«(0)| 2 ,(52) 



which is to be compared with p4|) . We see that 
both ratios (|13I17[) do not substantially depend on 
the mass of the radial excitation, and the measure- 
ment of the decay of i?-meson to the radially ex- 
cited D 1 should be feasible if the form factor f? (0) 
is significantly enhanced by the power corrections. 

To further emphasize the experimental feasibility 
of measuring B(B — > D'n) decays for large /+ (0), 
we can use /+(0) = 0.4 as obtained in ref. [l5| 
and compare this decay to B(B —> D^tt) that was 
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measured at Belle and BaBar in both cases, Class- 
I [H and Class-Ill 0. We get 



B{lf 



D'+w-) 



B(B° 
B(B~ 



> D'°TT- 



B(B- 



= 1.6(3) [1-2(3)] 



1.4(3) [1.1(3)], (53) 



where the first number on the r.h.s. corresponds 
to rriD' /mjj — 1.36 and the number within the 
brackets to mni/mo = 1.55(9). In other words 
the measurement of B(B — > D'n) is as feasible 
as that of B(B — > D^tt), provided the form fac- 
tor /+ (0) is large as recently suggested in the lit- 
erature. Note also that for the above estimates 
we used B(lf -> D; + tt~) = 4.9(7) x 1CT 4 , and 

B(lf -> D* 2 + tt-) = 8(1) x 10~ 4 Q. Experimen- 
tally, these measurements are made in the Dalitz 
plot analysis of B — > D*tttt decay. It is therefore 
important to mention that the statement made in 
eq. (|53|) remains valid even after we take into ac- 
count the results of refs. [3j, 152J, namely B(D' — > 
D*n) > B(D* -> D*tt). In other words, for a 
large value of the form factor /+ (0) the decay mode 
B — > D'tt should be as discernible in the Dalitz plot 
analysis of B -» D*irir as B — > D^tt. 



V. SUMMARY 

In this paper we develop the possibilities to check 
experimentally the size of the form factor f+(q 2 ) 
that parameterizes the B — > D' weak transition 
matrix element [5]. If its value is large, as recently 

claimed, the decay B — > D' + ir~ should be acces- 
sible experimentally and we predict its branching 
fraction in terms of |/+(0)| 2 , by using the mea- 
sured S(S° — > D + ir~) and by employing the fac- 
torization approximation to describe the decay am- 
plitudes. The uncertainties to the used approxima- 
tion are expected to be small since the ratios of the 
measured Class-I non-leptonic decays are known to 
be very well described by factorization. 

If the experimenters succeed in measuring 
B(B — > D' + n~), it would be interesting to check 
whether or not the corresponding |/+(0) | 2 is con- 
sistent with the ratio between the Class-Ill and 
Class-I B — > D'tt decay modes (fBTj) , which we could 
also predict using the factorization approximation. 
Note that the factorization in Class-Ill decays re- 
quires more assumptions as well as the computation 
of the decay constant fn 1 ■ 

We computed m£>'/m£) and Jd'/Id by using 
the gauge field configurations with Nf — 2 mass- 
degenerate light quark flavors, generated at four 



lattice spacings and for several light sea quark 
masses. We find, 



— = 1.55(9), tR. = 0.57(16). (54) 
niD Jd 



If the state observed by BaBar collaboration is 
indeed D', then our result is larger than theirs, 
rriD'/mD = 1-36. More research on both sides is 
needed to clarify the (potential) discrepancy. On 
the lattice QCD side it would be interesting to 
check whether or not tud 1 /tud becomes sensitive to 
the variation of the light quark mass in the region 
with very light quarks (closer to the chiral limit), 
the region not explored in the present study. Such 
a situation, that a hadron mass strongly depends 
on the sea quark mass when the latter is close to 
the chiral limit, was observed in the case of the 
Roper resonance on the lattice [4l[. Concerning 
the interpretation of the state observed by BaBar at 
2539(8) MeV, it is important to understand why its 
width is much larger than predicted. As a starting 
point one could verify if the predictions of ref. [3j] 
remain stable if one uses different set of wave func- 
tions (for example those of the model of ref. 0) or 
different models. 

We also improved the computation of the decay 
constants /d w by relying on the chiral and con- 
tinuum extrapolation of the ratios fD a /TriD, and 
(4>D s /4>D)/(fK/ 'fir)- More specifically we obtain: 



f Da = 252(3) MeV, 



/d. 



1.23(2)(1), (55) 



where the second error in the latter result re- 
flects the uncertainty due to inclusion/omission of 
the chiral logarithms in the light mass extrapola- 
tion to the physical limit. These two results give 
f D = 205(5) (2) MeV. 6 

Finally, our prediction for the non-leptonic de- 
cays depends only mildly on the choice of the mass 
of the radial excitation. For mjj'/fnD = 1-36 we 
have, 

B(B° -^D'+ir-) = [4.7(4) x 10- 3 ] |/f (0)| 2 , 
B(B- -^D'°7r-) = 



[5.1(5) x 10" 



0.13(4) 



l/£(0)| 2 , 
(56) 



6 Recent estimates of the vector meson decay constants /rj» 
Id* I 'Sd* i can be found in ref. |53|1 . 
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and for mw /mp = 1.55(9) we get, 



B(W -^D'+tt-) = [3.6(6) x 1(T 3 ] |/«(0)| 2 , 



[3.8(7) x 10" 



0.16(5) 



/f(0) 



I/4W 



(57) 
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p 


3.8 


3.9 


3.9 


4.05 


4.2 


4.2 


L s x T 


24 3 x 48 


24 3 x 48 


32 3 x 64 


32 3 x 64 


32 3 x 64 


48 3 x 96 


# meas. 


240 


240 


150 


150 


150 


100 


/^scal 


0.0080 


0.0040 


0.0030 


0.0030 


0.0065 


0.0020 


A 1 sca2 


0.0110 


0.0064 


0.0040 


0.0060 






/^sca3 








0080 






a [fm] 


0.098(3) 


0.085(3) 


0.085(3) 


0.067(2) 


0.054(1) 


0.054(1) 


Ma 


0.0194(7) 


0.0177(6) 


0.0177(6) 


0.0154(5) 


0.0129(5) 


0.0129(5) 




0.2331(82) 


0.2150(75) 


0.2150(75) 


0.1849(65) 


0.1566(55) 


0.1566(55) 



TABLE I: Lattice ensembles used in this work with the indicated number of gauge field configurations. Lattice spacing is 
fixed by using the Sommer parameter ro/a l5d . with tq = 0.440(12) fm fixed by matching /„■ obtained on the lattice with 
its physical value (c.f. ref. [32] )■ Quark mass parameters /i are given in lattice units. 



N{ 


P (csw) 


L 3 


x T 


# meas. 


^sca 




Kc 





6.2 (1.614) 


24 3 


x 48 


200 




0.1348 


0.125 


2 


5.4 (1.823) 


24 3 


x 48 


160 


0.13625 


0.1359 


0.126 



TABLE II: Lattice set-up for the results obtained by using the Wilson gauge and the Wilson-Clover quark action. K aea , k s 
and k c stand for the value of the hopping parameter of the sea, strange and the charm quark respectively. 
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(T ft a \ 


MS/o p vi 

llt'q y£ \JC V J 


m D q 


m D , q 


JD q 




(24, 3.80, 0.0080) 


0.0398(11) 


0.843(1) 


1.32(6) 


0.136(1) 


0.098(8) 


(24, 3.80, 0.0110) 


0.0547(15) 


0.852(1) 


1.34(3) 


0.139(1) 


0.105(6) 


(32, 3.90,0.0040) 


0.0216(5) 


0.741(1) 


1.05(5) 


0.110(1) 


0.054(8) 


(24, 3.90,0.0064) 


0.0345(8) 


0.7748(1) 


1.09(7) 


0.112(1) 


0.063(18) 


(24, 3.90,0.0085) 


0.0458(11) 


0.748(2) 


1.06(4) 


0.113(1) 


0.056(9) 


(24, 3.90,0.0100) 


0.0539(13) 


0.755(1) 


1.10(3) 


0.116(1) 


0.069(5) 


(32, 4.05,0.0030) 


0.0162(4) 


0.608(2) 


0.95(3) 


0.083(1) 


0.055(5) 


(32, 4.05,0.0060) 


0.0216(5) 


0.616(1) 


1.02(4) 


0.087(1) 


0.063(6) 


(32, 4.05,0.0080) 


0.0249(7) 


0.621(1) 


1.02(6) 


0.090(1) 


0.068(9) 


(32, 4.20,0.0065) 


0.049(2) 


0.521(1) 


0.79(3) 


0.071(1) 


0.038(4) 


(48, 4.20,0.0020) 


0.0150(7) 


0.497(1) 


0.81(5) 


0.064(1) 


0.037(7) 



TABLE III: Masses and decay constants, m D (i) and f D c), as computed from the solution to the GEVP discussed 
in the text. Note that the light valence quark and the sea quarks are degenerate in mass, m„ with the renormalized 
value given in the MS scheme. Note that the hadron masses and decay constants are given in lattice units while 
m^ s (2 GeV) is given in physical units [GeV], 
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IT ft a \ 

(Li, P, llq) 


TTlD s 


mD 's 


f r , 

JD S 


Jd' b 


f r, / f 

Jul Jit 


(24, 3.80, 0.0080) 


0.8703(8) 


1.36(7) 


0.1459(9) 


0.107(13) 


1.075(4) 


(24, 3.80, 0.0110) 


0.8717(7) 


1.35(3) 


0.1462(9) 


0.109(5) 


1.051(3) 


(32, 3.90,0.0040) 


0.7708(8) 


1.09(4) 


0.1206(7) 


0.062(7) 


1.114(7) 


(24, 3.90,0.0064) 


0.7715(8) 


1.11(4) 


0.1203(5) 


0.066(11) 


1.072(2) 


(24, 3.90,0.0085) 


0.7963(10) 


1.09(3) 


0.1199(7) 


0.064(8) 


1.057(2) 


(24, 3.90,0.0100) 


0.7713(8) 


1.12(2) 


0.1214(6) 


0.074(4) 


1.045(1) 


(32, 4.05,0.0030) 


0.6344(11) 


0.99(3) 


0.0923(7) 


0.063(5) 


1.129(6) 


(32, 4.05,0.0060) 


0.6355(9) 


1.03(4) 


0.0930(6) 


0.066(5) 


1.072(1) 


(32, 4.05,0.0080) 


0.6361(9) 


1.03(5) 


0.0941(8) 


0.072(8) 


1.050(2) 


(32, 4.20,0.0065) 


0.5243(7) 


0.81(2) 


0.0750(6) 


0.041(4) 


1.049(2) 


(48, 4.20,0.0020) 


0.5198(4) 


0.82(2) 


0.0726(3) 


0.042(3) 


1.134(6) 



TABLE IV: Similar as in Tab. I JJJ1 except that the valence quark mass is fixed to the strange quark mass value. We 
also list the values of Jk / fn obtained on each lattice (also referred to in the text as fx q /fw qq ) which are extracted 
in the same way as in ref. jT3] and corrected for the small finite volume effects f5l[ /. 



Quantity (Q) 


Fit form 




Bq 


Cq 


fD s /m,D s 




eq. 


p| 




0.1278(11) 


0.4(2) 


0.21(1) 


r — 4>D s /(t>D 


eq. 




[X = 


0] 


1.165(10) 


-0.34(3) 


0.006(5) 


r = 4>dJ4>d 


eq. 


cu 


[x = 


1] 


1.319(14) 


0.06(3) 


0.002(5) 


R = t/Uk/U) 


eq. 


Eg 


[x = 


0] 


1.018(6) 


-0.013(22) 


-0.002(4) 


R = t/Uk/U) 


eq. 


Eg 


[x = 


1] 


1.029(7) 


0.021(24) 


-0.002(4) 


m D i/m D 




eq. 


(|44l) 




1.55(9) 


0.6(1.5) 


-0.04(5) 


m D iJm Ds 




eq. 


(f44|) 




1.53(7) 


1.0(1.2) 


-0.05(4) 


fD>/.fD 




eq. 


(f44|) 




0.50(12) 


1.7(5.6) 


0.24(23) 


fD'Jfn. 




eq. 






0.52(9) 


1.5(4.6) 


0.19(21) 



TABLE V: Fit results of the quantities computed in this paper on the lattice. 



